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ABSTRACT 

We calculate the cooling times at constant density for halos with virial temperatures 
from 100 K to 1 x 10^ K that originate from a 3a fluctuation of a CDM power spectrum 
in three different cosmologies. Our intention is to determine the first objects that can 
cool to low temperatures, but not to follow their dynamical evolution. We identify 
two generations of halos: those with low virial temperatures, Tvir ^ 9000 K that remain 
largely neutral, and those with larger virial temperatures that become ionized. The 
lower-temperature, lower-mass halos are the first to cool to 75 percent of their virial 
temperature. The precise temperature and mass of the first objects are dependent 
upon the molecular hydrogen (H2) cooling function and the cosmological model. The 
higher-mass halos collapse later but, in this paradigm, cool much more efficiently once 
they have done so, first via electronic transitions and then via molecular cooling: 
in fact, a greater residual ionization once the halos cool below 9000 K results in an 
enhanced H2 production and hence a higher cooling rate at low temperatures than for 
the lower-mass halos, so that within our constant-density model it is the former that 
are the first to cool to really low temperatures. We discuss the possible significance 
of this result in the context of CDM models in which the shallow slope of the initial 
fluctuation spectrum on small scales leads to a wide range of halo masses (of differing 
overdensities) collapsing over a small redshift interval. This "crosstalk" is sufficiently 
important that both high- and low-mass halos collapse during the lifetimes of the 
massive stars which may be formed at these epochs. Further investigation is thus 
required to determine which generation of halos plays the dominant role in early 
structure formation. 

Key words: cosmology: theory — early Universe — galaxies: formation — molecular 
processes 



1 INTRODUCTION 

The study of the first generation of objects in the Universe 
that are able to cool sufficiently to collapse and form lumi- 
nous objects is a well-defined problem. By definition, there 
are no stars or other sources of ionizing radiation, and one 
does not have to consider feedback from supernovae and en- 
richment of the Universe with metals. The first objects to 
form arise from the collapse of high-sigma fiuctuations in 
the background density field. These peaks will virialize and 
begin to cool. Objects with virial temperatures T <, 9 000 K 
are cooled by H2 molecules: the molecules are excited by 
collisions with neutral hydrogen and then spontaneously de- 
excite with the emission of a photon. In the absence of 
metals the dominant cooling mechanisms for 9 000 K ^ T <, 
50 000 K are collisional excitation of neutral hydrogen and 
recombination of ionized hydrogen, and for T >, 50 000 K is 
collisional excitation of ionized helium. 
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The question of what are the masses of the first ob- 
jects to form in a standard CDM scenario was studied by 
Tegmark et al. (1997, hereafter T97) in a landmark paper 
entitled "How small were the first cosmological objects?" 
They analytically tracked a top-hat collapse to the point of 
virialization, at which point the gas was cooled at constant 
density. They accepted an object as having cooled if it met 
the criterion r(0.75zvir) < 0.75 Tvir, where Tvir is the virial 
temperature and z^ir the virialization redshift. They found 
that the first generation of objects that cooled in a stan- 
dard CDM scenario, virialized at a redshift of 27 and had 
a baryonic mass of about 10^ M©. In a later paper, Abel et 
al. (1998) redid the calculation with a difi'erent H2 cooling 
function and estimated a very similar virialization redshift 
but a smaller baryonic mass, 7 X 10^ Mq. 

In the present paper we adopt a similar approach to 
T97 but consider the collapse of not only small objects with 
virial temperatures <9 000K, but also objects of virial tem- 
peratures up to 100 000 K. The motivation for this is that 
the high virial temperatures will partially ionize the gas. 
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Since the gas will then cool more rapidly than it can recom- 
bine, the ionization level at temperatures <9 000K will be 
greater than it otherwise would have been if no re-ionization 
had taken place. This in turn accelerates production of H[2, 
ultimately resulting in enhanced cooling at lower tempera- 
tures. A similar effect was noted by MacLow & ShuU (1986) 
and Shapiro & Kang (1987) for the cooling of gas behind 
intergalactic shocks. 

We confirm the earlier results on the masses of the 
smallest halos that can cool, but show that these are not 
the first objects to do so in this constant-density model. 
A plot of cooling redshift versus halo mass shows two sepa- 
rate maxima corresponding to halos with virial temperatures 
of about 4 000K and 11 000 K. Thus there are two distinct 
generations of primordial halos depending upon whether or 
not their virial temperature is high enough to ionize the 
gas. Given our simplified model, each of these generations 
has a unique halo mass that will be able to cool first to 
low temperatures and form stars. If they virialized at the 
same time, then the higher-mass halos have the potential 
for more efficient cooling and, depending upon how the col- 
lapse occurs, could win out; however, in CDM cosmologies, 
the collapse of larger objects occurs slightly later than the 
collapse of smaller ones at the same overdensity relative to 
the rms at each scale. The enumeration of the difference in 
cooling times between the two generations of halos is the 
main topic of this paper. We show that this difi'erence is of 
order lOMyr, comparable to the lifetimes of massive stars. 
It is thus possible that the second generation of massive ha- 
los could collapse before feedback from the first stars affects 
their internal composition. The result of this paper, there- 
fore, is to argue that more detailed modelling is required to 
establish the relative importance of the two generations for 
subsequent structure formation. 

A second difference between our study and previous 
ones is that we consider more up-to-date cosmological mod- 
els. In particular, we choose a density fluctuation spectrum 
with less power on small scales. In a critical density universe, 
this has the effect of delaying collapse until a redshift of 11. 
Higher redshifts can be recovered by the introduction of a 
cosmological constant. 

We describe our chemical model in Section ^ and our 
numerical method in Section ^. The results are presented in 
Section ^ and discussed in Section Finally, we summarize 
our conclusions in Section ^. 



2 GAS CHEMISTRY AND COOLING 
FUNCTIONS 

In this section we present a minimal model of the gas chem- 
istry needed to accurately follow the temperature evolution 
of the gas. Our model is closely based on those of Abel et al. 
(1997) and Fuller & Couchman (2000) who both did a thor- 
ough search to identify the key reactions, but differs from 
theirs in the manner in which it handles the production and 
destruction of H~ and Hjj'. As a consequence, we end up 
with slightly different terms in our final equation for H2 
production. In principle this could be quite important; how- 
ever it seems to make little difference to the results over the 
parameter ranges considered in this paper. 

Table ^ lists the important reactions for the combina- 



tions of temperature and species abundances that we con- 
sider. The model is applicable to halos with temperatures 
between 100 K and 100 000 K and redshifts up to 50. Note 
that Hydrogen and Helium do not interact chemically. The 
latter is included in order that the ionization level be cor- 
rectly reproduced at high temperatures, but its omission 
would make very little difference to the results presented in 
this paper. We do not include photo-ionization from cosmic 
microwave background photons which is only important at 
redshifts in excess of 100. 



2.1 H2 production and destruction 

The equation for the rate of change of H2 abundance is de- 
rived in the appendix: 



dt 



R^R/i RqRy 
"o- + nu+ 



(1) 
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R9R8 



nuRio + n^-Rii 



where 

7? ^nuRr + n^-Rs, 



(2) 
(3) 



and R„ is the rate of Reaction n from Table 0. 

We next discuss the regimes under which each of these 
terms is important. 



2.1.1 Formation of H2 

The H~ channel is the dominant H2 formation path. The HJ 
channel is important only at high redshifts {z >, 200 when the 
H~ channel is suppressed by photo-destruction of H~that we 
have omitted here) or high temperatures. However, at high 
temperatures, the H2 is rapidly destroyed and for all the 
models that we consider in this paper, omitting the HJ for- 
mation channel makes no difference to the final abundance 
of H2. 



2.1.2 Destruction of H2 

It is the process of H2 destruction where the present work 
differs slightly from that of Abel et al. (1997) and Fuller & 
Couchman (2000). In these two papers, the term for destruc- 
tion by H+ is given as 



nH2"-H+-R9i 

whereas we have 



"-H2"'H+-R9- 



-Rs 



Rtj 



■ nH2?lH+-R9 



Rtj + 



(4) 



(5) 



This takes into account that, at low ionization levels, much 
of the HJ that is produced will be immediately converted 
back to H2 so that there will be no net destruction. In some 
of the models that we consider, this can make transitory 
differences of a factor of 10 or more in the H2 abundance. 
However, as H2 production takes over from destruction at 
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Table 1. This table summarizes the important reactions needed in order to calculate accurately the abundance of H2. References are: 
HTL Haiman, Thoul & Loeb (1996); GP Galli & Palla (1998); FC Fuller & Couchman (2000); SLD Stancil, Lepp & Dalgarno (1998); 
AAZN Abel, Anninos, Zhang & Norman (1997). 



Reaction 



Rate / cm^ s 



Reference 
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He+ + e- 
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2H 
-H++H 
■ 3H 
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He++ + 2e- 
I — » He+ + 7 



5.9 X 10-"rO-5(l + rO-5)-iexp(-1.58/T5) 

3.3 X 10-10 T-0-^(1 + T0-7)-i 

1.4 X 10-i**T0-93exp(-T4/1.62) 

1.3 X 10~9 

4.0 X 10-*^ T~°-^ 

2.1 X 10-23 ri-8exp(-2/Ti) 

6.4 X 10-^° 

1.2 X 10-^ Tg-o-* 

min(3.0 X IQ-^O cxp(-2.11/T4), 1.5 X 10-i°cxp(-l, 
7.1 X 10-"r2-0i(l + 2.13T5)-3-5iexp(-5.18/T4) 
4.4 X 10-i°r0-35exp(-1.02/r5) 
2.4 X 10-" (1 + rl"-^)-! exp(-2.85/r5) 



4O/T4)) 



1.5 X 10" 



+ 1.9 X 10" 



' exp(-5.64/T5) (0.3 + exp(9.40/T4)) 



5.7 X 10-12 (1 + TgO-S)-! exp(-6.32/r5) 
1.3 X 10-ST-°-^(1 + T0-7)-i 
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low temperatures, the effect on the final H2 abundance is 
quite small (at most a few percent). 

With the correct rate for the destruction channel, 
there is a small parameter range for which HIO becomes the 
dominant destruction process. However, this is so fleeting 
that it makes a negligible difference to the results and can 
safely be omitted. 

We note that the previous papers consider only low tem- 
peratures, T < 6 000K for which H2 destruction is relatively 
unimportant compared to its production, and the error in 
their results is negligible. Nevertheless, a minimal model 
that includes H2 destruction should use our Equation ^ 



2.2 Ionization level 

In principle, the Equation for the rate of change of abun- 
dance is every bit as complicated as that for H2. However, 
it turns out in practice that there are only two important 
terms: 



dt 



Similarly, for the Helium species 
driHc 



dt 

dnHc++ 



dt 



(6) 



(7) 



(8) 



The number density of electrons is given by charge conser- 
vation: 



nH+ + nHe+ + 27lHe + + • 



(9) 



al. (1996) and Forrey et al. (1997) which together cover a 
wide temperature range: 



lof 



AH,(r) 



1" ' nnnHjergcm^ s" 



-103.0 + 97.59 Tiog 



48.05 + 10.80 Tw 



(10) 



[log -f- lu.ou Jlog — 0.9032 Tiog, 

where Tiog = logio(r/K). 

For all other cooling processes, we use the rates given 
in Haiman, Thoul & Loeb (1996). The most important of 
these are as follows. At temperatures, T>, 10 000 if, coUi- 
sional excitation and (less importantly) ionization of atomic 
Hydrogen take over from molecular Hydrogen cooling: 

Ah,cc(T) 



Wg-nnerg cm'^ s~ 



AH.ci(T) 
riHcrgcm^ s" 



= 7.50 X 10" 



= 4.02 X 10" 



_Tl_ 

l + d 



(11) 



(12) 



where Tn is the temperature in units of 10" K. 

At even higher temperatures, 50 000 A" < T < 100 000, 
coUisional excitation of He"*" is the dominant coolant. 

-0.397 



A 



Hc+, 



17 ^0 



4 737 

5.54 X 10"''^^^ . (13) 

l+T/ 

We include other processes from Haiman, Thoul & Loeb 
(1996) in our model, including inverse Compton cooling from 
cosmic microwave background photons, but at the modest 
redshifts (z < 50) and high densities of the collapsed halos 
that we consider in this paper none contribute at more than 
the few percent level. 



^e- "■Hc+ erg cm'^ s^ 



2.3 Cooling Terms 

At low temperatures the main coolant is molecular hydro- 
gen. We use the cooling rate given in Galli & Palla (1998) 
as summarised by Fuller & Couchman (2000). This gives a 
fit to the low-density limit of the calculations of Martin et 



3 NUMERICAL PROCEDURE 
3.1 General approach 

The approach that we have adopted for this paper is to 
specify an initial temperature, ionization level and H2 frac- 
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tion and to follow their evolution at constant density, using 
the relevant cooling terms and reaction rates discussed in 
the previous section. The assumption of constant density is 
only valid for halos whose cooling times are much longer 
than their dynamical times. It is adequate for the purposes 
of this paper to assess if a given halo at constant density 
can begin to cool efficiently, but not to follow the collapse 
of halos once they do cool. 

The equations are integrated using the RK4 integrator 
from Press et al. (1992) modified to use an adaptive timestep 
that allowed neither the abundances nor the temperature to 
vary by more than 0.1 per cent during a timestep. We also 
tried a Bulirsch-Stoer integrator which gave identical results 
for a slightly poorer performance. For tests, we integrated 
simplified networks of equations for which a solution can be 
obtained analytically: for example. Equations 15 and 16 of 
T97. 

The initial halo parameters are picked to represent the 
conditions of a virialized object which has collapsed from a 
3(T peak in a CDM scenario, with virial temperatures, Tvir 
in the range 100 K to 100 000 K. We consider two different 
measures of cooling. Firstly, we measure how long it takes 
halos to cool to To. 75 = 0.75 Tvir — T97 define an object to 
have cooled if its temperature decreases by 25 per cent or 
more in the time that redshift does likewise. We find two 
distinct populations of clouds: low-mass ones that cool via 
molecular hydrogen and high-mass ones that cool via elec- 
tronic transitions. Secondly, we look at the amount by which 
an object can cool in one dynamical time. Clouds with virial 
temperatures greater than 9 000K (i.e. those that have re- 
ionized) form later (at the same o) but cool more effectively 
than lower-mass clouds, and have a much greater H2 fraction 
once they have cooled to low temperatures. 

3.2 Cosmological Models 

We compute the collapse redshift of objects arising from 3a 
peaks of a CDM power spectrum. The spectrum was calcu- 
lated using a real-space top-hat window function, the trans- 
fer function of Bond & Efstathiou (1984) for scales above 
1 /i~^Mpc and the transfer function of Bardeen et al. (1986), 
or BBKS, for smaller scales. We choose this combination as 
the Bond & Efstathiou transfer function is more accurate 
than BBKS but makes no attempt to accurately calculate 
the function on scales below 1 /i~^Mpc. 

We present results for three different cosmological mod- 
els, as listed in Table For comparision with previous work, 
we first choose the standard, cold dark matter cosmology 
(SCDM) with power spectrum shape parameter, F = 0.43. 
When normalised to the COBE results, this power spectrum 
is now known to have too much power on small scales and so 
the other two models that we consider use a smaller value, 
r = 0.21. 

The ACDM model is our best-guess at the most 
favourable CDM model. The combination of fio, and h 
naturally produces a power spectrum with the correct value 
of r, and the Hubble parameter and baryon fraction both lie 
close to currently preferred values (see, for example, Freed- 
man et al. 2001 for the former and Ettori & Fabian 1999 for 
the latter). The normalisation of the power spectrum, ag, 
is chosen to reproduce the correct local abundance of rich 
clusters (Viana & Liddle 1996). 



The rCDM model is an attempt to salvage a critical 
density model in the post-COBE epoch. It has some motiva- 
tion in decaying neutrino scenarios (e.g. Hansen & Villante 
2000) but is mostly phenomenological in nature. For this rea- 
son, we have fixed F to have the same value as in the ACDM 
model, but have used a lower value of ag to again reproduce 
the rich cluster abundance. In addition, we have had to take 
a lower value of the Hubble parameter in order to overcome 
the age problem, and a baryon density considerably higher 
than that predicted by primordial nucleosynthesis to give 
the correct baryon fraction in clusters. 

3.3 Initial halo properties 

We calculate the initial properties of halos by assuming that 
they have settled down into virial equilibrium. This will be 
a good approximation for halos whose cooling time substan- 
tially exceeds their dynamical time and should therefore en- 
able us to identify objects that are able to cool with rea- 
sonable accuracy. To correctly model the properties of ha- 
los whose cooling time is shorter than their dynamical time 
would require a more sophisticated model. 

In an isothermal sphere, the density is strongly peaked 
towards the centre of the halo. Here we follow the evolution 
of gas which has a density fixed at the mean value within 
the virialised halo, 

Pvir= (^)p60(l + 2:vir)^ (14) 

where Zvir is the virialization redshift; pi,o is the current 
mean density of baryons in the Universe; and Ac is the 
mean overdensity relative to the critical density within the 
virialized halo which we take to be 187r^ for the SCDM and 
rCDM cosmologies, and IStv'^QP-'^^ for the ACDM cosmol- 
ogy, where Q. is the density parameter at the time of virial- 
ization (Eke, Navarro & Frenk, 1998). We take a hydrogen 
mass fraction (in all its forms: neutral, ionized and molecu- 
lar) oiX = 0.755. Thus 

riH.tot = nn + + 2nH2 = pvirX. (15) 
Similarly, 

nHc,tot = riYlc + n-jic+ + "■Hc++ = PvirY, (16) 

where Y = {1 — X) = 0.245 is the Helium mass fraction. 

For the rCDM and ACDM cosmologies, the residual 
ionization level from the early Universe is taken from Peebles 
(1968, 1993) but divided by two as an attempt to compen- 
sate for the neglect of stimulated recombination as suggested 
by the results of Grachev & Dubrovich (1991) and Sasaki & 
Takahara (1993). Although dividing by two is somewhat ar- 
bitrary, it does result in a more accurate estimate of the 
initial ionization level than many previous papers without 
explicitly solving the time-dependent evolution of the ion- 
ization level with the inclusion of stimulated recombination. 
We take 

^^ = 6xlO-«fi|(nM,ft)-\ (17) 
'^H,tot 

which gives the residual ionization levels listed in Table ^ 

For halos hotter than about 8 000-9 000 K, the gas 
will be ionized above the residual value by the pro- 
cess of virialization. For these we use the equilibrium 
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Table 2. Cosmological parameters for the three models: model name; density parameter; cosmological constant in units of Aq = A/3H^; 
Hubble parameter in units of h = Hq/ 100km s~^Mpc~^; root-mean-square dispersion of the density within spheres of radius 8 /i~^Mpc; 
residual ionization fraction. 



Name 


no 


A(, 




h 


r 




ne,ros/nH,tot 


SCDM 


1.0 


0.0 


0.076 


0.5 


0.43 


0.60 


1.58 X lO--* 


rCDM 


1.0 


0.0 


0.184 


0.5 


0.21 


0.60 


6.52 X 10-5 


ACDM 


0.35 


0.65 


0.038 


0.7 


0.21 


0.90 


1.33 X 10"* 



values determined by solving Equations ^ ^ and ^ for 
dnjj+ /dt =dnuc/dt =dnHe++ /dt = 0: 



Ri 



R1 + R2 



R12R15 



nn.tot 

n-Hc,tot 

and 

Whc++ _ 

?^Hc,tot -R13-R15 + R12R14 + R12R15 



R13R15 + R12R14 + i?12i?15 



R12R1 



(18) 
(19) 

(20) 



The background level of H2 in the Universe is taken to 
be /h2,„, = 1.1 X lO"*^ as calculated by Galh & Palla (1998). 
We take the initial temperature to be 

finiH GMtot 



40.8 ■ 



2rvi, 



1.225 



1 + Zvir 



Aft, 



105 Mo 



21) 



K. 



Here Mtot is the total mass (dark plus baryonic) which we 
assume is distributed as an isothermal sphere within the 
virial radius, r^ir; mn is the mass of a hydrogen atom; ks 
is the Boltzmann constant; G is the gravitational constant; 
and /X is the mean mass of particles in units of mn'- 

■n-H.tot + 4nHc,tot 



M = 



nH,tot — + n-Hc 



(22) 



When the molecular Hydrogen abundance is low, then /i 
varies between 1.23 at low temperatures, r;S10 000K, and 
0.59 at very high temperatures, T > lOOOOOK. 

We shall use the halo virialization redshift, .Svir , as Our 
ordinate in most of the plots that follow. This can be con- 
verted to temperature or mass in a cosmology-dependent 
way. We plot these relations in Figure |l] and list examples of 
the conversion between temperature and other quantities in 
Table ^ The kink in the temperature profiles corresponds 
to the changing value of fi from Equation ^ during Hydro- 
gen ionization (there is a second kink due to ionization of 
Helium but that is barely discernible). Note that halos of a 
given temperature have similar total masses in each cosmol- 
ogy, although they may virialize at very different redshifts. 



4 RESULTS 
4.1 SCDM 

To compare with previous work, we first consider the col- 
lapse of halos in the SCDM cosmology. We do this for two 
different forms of the H2 cooling function. The first, de- 
scribed above in Equation 0, is from Galli & Palla (1998, 
hereafter GP98); the second is from Lepp & ShuU (1984; 



Table 3. Conversion between virialization temperature, redshift 
and mass for the three models: temperature; model; virialization 
redshift; total mass; baryonic mass. 



T/K 


Model 




A/tot/ M0 


A/bary/Mg 


100 


SCDM 


31.6 


4.2 


X 


103 


3.2 X 10^ 




tCDM 


18.0 


9.2 


X 


10^ 


1.7 X 10^ 




ACDM 


34.0 


2.7 


X 


10^ 


2.9 X 102 


1000 


SCDM 


25.0 


1.8 


X 


10'> 


1.4 X lO"* 




rCDM 


14.3 


4.1 


X 


10'> 


7.5 X lO-* 




ACDM 


27.3 


1.1 


X 


10^ 


1.2 X 10* 


10 000 


SCDM 


18.9 


8.7 


X 


10" 


6.6 X 10^ 




rCDM 


10.9 


1.9 


X 


10^ 


3.5 X 10*^ 




ACDM 


21.0 


5.4 


X 


lO'' 


5.9 X 10^ 


00 000 


SCDM 


11.6 


1.6 


X 


109 


1.2 X 10^ 




rCDM 


6.8 


3.3 


X 


109 


6.1 X 10^ 




ACDM 


13.3 


9.7 


X 


108 


1.1 X 10^ 



hereafter LS84). This latter form is the one used by Abel 
et al. (1998); it gives an order of magnitude more cooling 
at temperatures below 1000 K and so favours the collapse of 
low-mass, low-virial temperature objects. 

Initially the gas within the halos is assumed to be 
shocked to the virial temperature and to be pressure sup- 
ported. If it can cool significantly within one dynamical time 
then it will contract towards the centre of the halo in order 
to maintain pressure support. Accordingly, we plot in Fig- 
ure ^ the redshift, 2:0.75, at which the gas (if it maintains 
constant density) will lose 25 percent of its initial energy. 
The solid curve was generated using the GP98 cooling func- 
tion and the dashed curve using that of LS84. The lower 
dotted line shows the relation 20.75 = 0.75 Zvir which is the 
condition used by T97 to separate clouds which can cool 
from those which cannot. The dashed line shows the redshift 
one dynamical time, tdyn, after the time of virialization, t^ir, 



where we set tdy 



First concentrate on the GP98 curve. This intersects 
the line 20.75 = 0.75 Zvir at 2vir ~ 23.8 which means that the 
smallest clouds that can cool according to the T97 criterion 
will have virial temperatures of rvir=l 600 K and masses of 
4 X 10^ Mq . The most important quantity for comparison 
with the results of T97 is the temperature which agrees very 
well with that given in their Figure 5. For a virialization 
redshift of 23.8, they require a slightly higher temperature 
of 2 500 K for collapse, but this difference is attributable to 
our having a slightly higher baryon density and H2 cooling 
rate. 

T97 assume, within the constraints of the constant den- 
sity model, that the smallest halos that can cool will be the 
first to do so, but that is not the case. As one moves to 
lower 2vir and higher Tvir, the efficiency of H2 production 
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Figure 1. This figure sliows the relation between virializa- 
tion redsliift and (a) temperature and (b) total mass for the 
three cosmologies: SCDM (solid line), rCDM (dashed line), and 
ACDM(dotted Une). 



increases and the cooling time decreases. Consequently, the 
first objects to cool to 75 per cent of their initial temperature 
virialise later at Zvir ~ 21.6, and have higher virial temper- 
atures (3 700K) and masses (1.6 x 10^ M©). This maximum 
in 20.75 corresponds roughly to the redshift one dynamical 
time after virialisation. 

At higher temperatures the efficiency of H2 cooling con- 
tinues to increase but cannot compensate for the later viri- 
alization redshift and so zo.75 decreases once more. Note 
the cusp in the curve at Zvir ~ 19.2, corresponding to 
Tvir ~ 9 000 K. As this temperature is approached, H2 starts 
to be destroyed and so its abundance plummets and the 
cooling time increases. However, at higher temperatures col- 
lisional excitation and ionization take over and these are 
much more efficient cooling processes so the cooling time 
drops once more. In fact for temperatures above 11 000 K 
the cooling is essentially instantaneous (i.e. the cooling time 
much less than the dynamical time) . 

Thus the two peaks in Figure H correspond to two dif- 



Figure 2. Redshift at which a gas halo will cool to 75 per cent 
of its initial temperature, 2o.75i versus virialization redshift, 2viri 
for cosmology SCDM for two different cooling functions: dashed 
curve, LS84; solid curve, GP98. The upper and lower dotted lines 
show the relations 2:0.75 = -Zvir and zo.75 = 0.75 Zvir, respectively; 
the dash-dotted line shows the redshift one dynamical time after 
virialization. 



ferent classes of object: the first to collapse, at ^0.75 = 19.4, 
are dominated by molecular cooling and the second, at 
Zo.75 = 18.5 by electronic transitions. As we are interested in 
the first objects to cool, it might be thought that the smaller 
halos are the more important, but this is not entirely clear 
because the T97 criterion does not trace the cooling down to 
very low temperatures. In fact the more massive halos have 
a large residual ionization which results in a greater H2 pro- 
duction once they have cooled below 9 000K, and so they 
may be the first objects to cool to rea lly l ow temperatures. 



below. 



We will discuss this further in Section 4.2, 

Turning now to the LS84 curve, we see that this is qual- 
itatively similar. However, it predicts shorter cooling times 
for low-mass clouds because the H2 cooling rate is much 
higher. Hence the smallest objects that can collapse have a 
lower mass and temperature than for the GP98 cooling func- 
tion. The lowest-mass objects that satisfy the T97 criterion 
virialize at a redshift of 26.6 and have virial temperatures 
of 570 K — this seems to agree reasonably well with the pre- 
dictions from Figure 12 of Abel et al. (1998). The peak of 
the curve has moved to Zvir ~ 24.7 and the temperature and 
mass of these first objects are significantly lower, 1 100 K 
and 2.3 x 10* Mq, respectively. By contrast, the properties 
of the halos corresponding to the higher-mass peak are only 
slightly modified. 

We use the more up-to-date GP98 cooling function 
throughout the rest of this paper. 



4.2 rCDM 

We now turn to the rCDM cosmology. Like the SCDM cos- 
mology, this has a critical density of matter and a high 
baryon density, but it has a spectral shape that gives less 
power on small scales and hence the first objects collapse 
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Figure 3. The solid line shows the redshift at which a gas halo 
will cool to 75 per cent of its initial temperature, ^o.75 > versus viri- 
alization redshift, Zvir, for cosmology tCDM using the GP98 cool- 
ing function. The upper and lower dotted lines show the relations 
•20.75 = 2vir and 20.75 = 0.75 Zvin respectively; the dash-dotted 
line shows the redshift one dynamical time after virialization. 



Figure 4. Temperature obtained after one dynamical time versus 
virialization redshift, for cosmology tCDM. The dotted line shows 
the initial, virial temperature of the halo. The cusp in the curve 
corresponds to a virial temperature of 9 000 K. 



at much lower redshift. This is shown in Figure ^ which is 
the analogue of Figure ^ but for the rCDM cosmology. We 
can see that gas halos first manage to cool to 0.75 Tvir only 
at redshift «o.75 ~ 10.8: these have virial temperatures of 
4 500K and total masses of 5.0 x 10'=' M©. 

In order to form stars, the gas has to cool to much 
less than To. 75. Therefore we consider a second measure of 
cooling: Figure ^ shows the final temperature of the halo 
after one dynamical time (once again assuming constant 
density — this will be a good assumption only for clouds 
which cool by a small amount in a dynamical time; for oth- 
ers it will underestimate their cooling rate). The coldest 
clouds are those to the right-hand-side of the Figure, but 
this is only because they were born with low temperatures: 
they have cooled very little. The final temperature is an in- 
creasing function of virial temperature up to about 9 000 K 
when it suddenly plummets, so that somewhat surprisingly, 
the larger, higher virial temperature halos, have cooled to 
a lower temperature than the smaller ones. The reason for 
this is shown in Figure ^ which shows that the H2 abun- 
dance after one dynamical time is larger in high-mass halos 
than low-mass ones. This is because they become highly ion- 
ized and their residual ionization once they have cooled back 
down below 9 000K is greater than that of low-mass halos, 
which in turn results in a greater production of H2. Note 
that the cooling rate of these high-mass halos will continue 
to be greater than that of low-mass halos even after one dy- 
namical time. It is therefore conceivable that they will be 
the first to cool to really low temperatures. Because the gas 
is free to move around on this timescale, however, realistic 
dynamical simulations are required to model the problem 
accurately. 




Figure 5. Ionization fraction (solid line) and H2 abundance 
(dashed line) after one dynamical time, for cosmology rCDM. 
The dotted line shows the initial ionization level of the halo. 



4.3 ACDM 

This cosmology is currently the most popular CDM model. 
Halos virialize at higher redshift than in rCDM but qualita- 
tively, their behaviour is just the same. Figures ^-|^ mimic 
those of Figures ^|-^ and show all the same features. In this 
cosmology the first halos to cool to To. 75 have virial temper- 
atures of 3 400K, masses of 8.9 x 10^ Meand cool to 75 per 
cent of their virial temperature at 20.75 ~ 21.8; the second 
generation of higher mass halos cool at z ~ 20.6 and have 
virial temperatures and masses of 10 600 K and 5.9 x 10® Mq . 
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Figure 6. The solid line shows the redshift at which a gas halo 
will cool to 75 per cent of its initial temperature, zo.75i versus viri- 
alization redshift, Zvir, for cosmology ACDM using the GP98 cool- 
ing function. The upper and lower dotted lines show the relations 
■20.75 = -Zvir S'lid ^0.75 = 0.75 2:vir> respectively; the dash-dotted 
Una shows the redshift one dynamical time after virialization. 




Figure 7. Temperature obtained after one dynamical time ver- 
sus virialization redshift, for cosmology ACDM. The dotted line 
shows the initial, virial temperature of the halo. 



5 DISCUSSION 

The properties of the halos corresponding to the two peaks 
in Figures ^, ^ and |^ are summarized in Table ^. 

Note that, although the redshift at which the first ob- 
jects cool is quite different in different cosmologies, their 
virial temperatures are very similar, reflecting the similar 
chemical processes that are occurring within them. 



Figure 8. Ionization fraction (solid line) and H2 abundance 
(dashed line) after one dynamical time, for cosmology ACDM. 
The dotted line shows the initial ionization level of the halo. 



Table 4. Properties of the two generations of halos to form in 
each of the cosmologies (corresponding to the molecular hydrogen 
and electron-cooling peaks of the curves in Figures ^, ^ and^: 
cosmological model; redshift at which the halo cools to 75 per 
cent of the virial temperature; virial temperature; total mass of 
the halo; baryonic mass of the halo. 



Model 


2^0.75 


Tvir/K 


Mtot/Mg 


A4ary/Mc 


Generation 1 










SCDM 


19.4 


3 700 


1.6 X lO** 


1.2 X 10^ 


rCDM 


10.8 


4 500 


5.0 X 10^ 


9.1 X 10^ 


ACDM 


21.8 


3 400 


8.9 X 10^ 


9.7 X 10* 


Generation 2 










SCDM 


18.5 


10 800 


9.9 X 10^ 


7.6 X 10^ 


rCDM 


10.7 


10 600 


2.1 X lO'^ 


3.8 X 10^ 


ACDM 


20.6 


10 600 


5.9 X lO** 


6.4 X 10^ 



5.1 The effect of cosmology 

Although it is often not stated explicitly, previous work 
seems to concentrate on the SCDM cosmology which is now 
known to be a poor model of the Universe; the currently 
favoured model being ACDM. Coincidentally these two give 
similar formation redshifts for the two generations of halos 
and have similar baryonic masses. The phenomenologically- 
based critical-density rCDM cosmology, on the other hand, 
forms its first objects at much lower redshifts, but never- 
theless still comfortably above the observed redshifts of the 
most distant objects yet observed. 

The baryonic masses of these first cooled halos are a 
factor of 10 times larger in the rCDM than in the ACDM 
cosmology — in the former they correspond roughly to the 
mass of globular clusters. Note, however, that the peaks cor- 
responding to the first generation halos in Figures |2|, ^ and 
^ are very broad so that we would expect halos with a wide 
range of masses to cool almost simultaneously. 
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5.2 The H2 mass fraction required for collapse 

We have shown in Section that the virial temperature 
required for coUapse agrees with the results from T97 and 
Abel et al. (1998). These two papers both claim that the H2 
fraction required for collapse of the first halos is 5 x 10~^, 
a result confirmed by Fuller & Couchman (2000). However, 
the H2 fractions shown in Figures gandg are much lower, 
about 1 X lO""*. There is no contradiction here, however. 
The Figures show the H2 fraction after one dynamical time, 
whereas T97 give the asymptotic H2 fraction at late times — 
for which we get a similar value of 4.4 x 10~^ for the SCDM 
cosmology (similarly 4.2 x 10~* for ACDM and a slightly 
lower value of 3.2 x lO""* for rCDM). 

5.3 The relative roles of the two generations of 
halos 

When talking about the first cosmological objects, we mean 
all those halos that can collapse and form stars before pol- 
luted by feedback of metals (and energy) from supernovae. 
Table ^ lists the differences in collapse times for the two 
generations of halos in each cosmology; these are similar to 
the life-time of the most massive stars (a 15 solar mass star 
has a lifetime of 10 Myr) . 

The time dilferences in Table ^ are likely to be upper 
limits because the enhanced H2 production will lead to more 
rapid cooling at low temperatures for the second generation 
halos. The question arises as to whether, even if small ob- 
jects are the first to collapse, larger objects may collapse 
around them and form stars before the first generation of 
supernovae explode. 

The results presented in this paper are for 3 a den- 
sity fluctuations. In reality there will be fluctuations with 
a continuous range of overdensities, leading to an enormous 
spread in collapse times. Note also that for the masses we are 
dealing with, the power spectrum is almost flat — in this sit- 
uation Press-Schecter theory (Press & Schechter 1974) pre- 
dicts that there is as much collapsed mass in each decade of 
mass as any other. 

We have taken a simple model in which halos have no 
substructure. In reality, high-mass object will contain larger 
fluctuations on smaller scales. For example, the inclusion of 
density profiles may raise the core density of even small ha- 
los to a point where free-fall collapse could proceed. In this 
case, we should also consider the effect of ionizing radiation 
from high-mass stars; however the amount of ionizing radi- 
ation and the degree to which this permeates the halo will 
be highly-dependent upon the initial-mass-function and the 
physical conditions in the star-forming region. 

The halo-in-halo problem can be addressed using a sta- 
tistical model of the collapse of halos but a realistic estimate 
of the star-formation time and the structure of star-forming 
regions requires numerical simulations of the dynamics of 
the collapse as discussed in the next section. 

5.4 Towards more realistic models 

Our model does not attempt to follow the increasing den- 
sity in objects that do cool on less than a dynamical time. 
To do so requires the use of some form of hydrodynamics 
code. Several groups have attempted this. Haiman, Thoul 



Table 5. The redshift, zo.75 and age, to. 75, of the universe at the 
time at which the first halos in each generation have cooled to 
75 per cent of the virial temperature. The final column gives the 
difference in age between the two generations. 

Model 20.75 to.75/Myr 20.75 to.75/Myr At/Myr 



Generation 1 Generation 2 



SCDM 


19.4 


140 


18.5 


151 


10 


tCDM 


10.8 


322 


10.7 


326 


4 


ACDM 


21.8 


145 


20.6 


157 


12 



& Loeb (1996), Omukai & Nishi (1998) and Nakamura & 
Umemura (1999) have all simulated spherically-symmetric 
collapses; Fuller & Couchman (2000) go further and simu- 
late a top-hat collapse with and without substructure using 
a three-dimensional code. In general, these results confirm 
the simple analytic predictions. 

One of the most important simplifications of our model 
is that it ignores substructure. Any halo that can cool and 
collapse on a dynamical time will contain smaller fiuctu- 
ations of even higher overdensity. This will result in frag- 
mentation of the cloud as the collapse proceeds. This has 
been investigated by Abel et al. (1998) and Abel, Bryan 
& Norman (2000) who have used a hierarchical grid code 
to perform very high-resolution simulations of the collapse 
and fragmentation of the first objects, down to scales of just 
1 pc. They find a very filamentary structure develops with 
the first stars forming from small knots at the intersection 
of filaments. Only a small fraction of the gas reaches suffi- 
ciently high densities to allow star-formation. If these first 
stars were to feed back energy into the surrounding cloud 
and disrupt it, then that would suggest that the size of the 
first star clusters may be much smaller than the size of the 
cloud from which they formed. However, the masses of these 
star clusters in the simulations is small and may not pro- 
duce many high-mass stars. Also, the simulations have not 
yet been run for long enough to see whether the surrounding 
cloud will collapse before feedback becomes effective, as we 
have suggested in Section |5.3| . It should be noted that very 
different results were obtained by Bromm, Coppi & Larson 
(1999) using a particle-based hydrodynamics, but with much 
poorer resolution. They found that the cloud collapsed to a 
rotationally-supported disk which then broke up into very 
massive star clusters. It will be a while, we suspect, before 
any concensus emerges. 

Fuller & Couchman (2000) simulated a cubical region 
of side 25 /i~^kpc in the SCDM cosmology with several dif- 
ferent realisations of the density fiuctuation spectrum, using 
an N-body, hydrodynamics code. This random realisation is 
just what is required to look at the relative importance of 
halos of different mass. They found that the most massive 
objects that collapsed within the region did so at a wide 
range of redshifts, 15-30, and had a similar large spread in 
mass. Unfortunately, their simulation volume was not large 
enough to sample our second-generation halos. 



6 CONCLUSION 

In this paper we have considered the cooling of gas within 
spherical, virialized halos in the high-redshift Universe. Our 
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technique is similar to that used by T97, but with a more 
up-date cooHng function and cosmological modeL In addi- 
tion, we have investigated halos with a wider range of virial 
temperatures, in the range of 100 K to 100 000 K. 

We have followed the abundances of Hydrogen and He- 
lium species, including molecular Hydrogen. We use a more 
complete destruction term for H2, taking into account the 
relative importance of the reverse reaction H7 which other 
authors have ignored. This makes a difference of a few per 
cent to the final H2 abundances for halos which cool from 
above 9 000K. 

The main coolant for temperatures below 9 000 K is 
molecular hydrogen. Unfortunately, the cooling rate seems 
poorly-known and a recent determination by Galli and Palla 
(1998) gives values at temperatures below 1 000 K over an or- 
der of magnitude lower than those of Lepp and ShuU (1984). 
Consequently the first objects to cool in the former do so 
later and have much higher masses and virial temperatures. 

We follow T97 by defining clouds to have cooled if they 
lose 25 per cent or more of their energy in the time that 
the redshift has decreased by 25 per cent. We obtain similar 
virial temperatures for the smallest halos that can cool, but 
note that the first objects to cool in each cosmology are 
more massive and have higher virial temperatures of about 
4 000 K. In the SCDM and ACDM cosmologies the formation 
redshift is 20.75 ~20; for the rCDM cosmology it is much 
lower, Z0.75 ~11- 

We identify a second generation of halos that cool about 
10 Myr after the first one. These are halos with virial tem- 
peratures in excess of 9 000 K for which there is a signifi- 
cant fraction of free electrons. The cooling is dominated by 
electronic transitions at high temperatures and is almost in- 
stantaneous, occurring on much less than a dynamical time. 
Just as significant, however, is the fact that the residual 
ionization is greater than in low-mass halos and the produc- 
tion of H2 is much greater. Consequently, they cool to much 
lower temperatures in a dynamical time than do the first 
generation of halos. 

Our model suffers from three deficiencies: we consider 
only 3(7 density fiuctuations, we ignore substructure, and 
we cannot follow the collapse of halos whose cooling times 
are shorter than their dynamical times. More sophisticated 
studies are required to determine whether the first or second 
generation halos are more important for determining the 
mass of the first star-clusters, or whether they both have a 
role to play. We hope to report on these in future papers. 
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The appendices describes two different ways of deriving 
Equation hi for the rate of change of abundance of H2. We 
do this because nowhere in the literature does it seem to be 
spelt out in detail and because our equation differs slig htly 
from those used previously, as described in Section |2.l[ 



APPENDIX A: DERIVATION VIA REMOVAL 
OF H AND H+ FROM THE REACTION 
NETWORK 

The reaction rates for destruction of H~ and Hj are much 
greater than for formation or destruction of H, H^, e~ and 
H2. Hence we can regard H~ and Hj as short-lived species 
that are destroyed the instant that they are produced. This 
has the advantage that we can eliminate them from the re- 
action network, as described below. 

Consider first H~. This has two destruction channels, 
Reactions H4 and H5 in Table |l|, whose combined rate is 

R^- =nu_RA + n^+Rrj, (Al) 

where Ra and i?5 are the reaction rates as listed in the Table 
(at redshifts greater than 110 photo-destruction of H~ by 
cosmic microwave background photons is also important, 
but we do not consider such high redshifts in this paper). 
Hence the fractions of H^ that decay via reactions H4 and 
H5 are nuRt/ Rh- a-nd njj+ R5/R11- , respectively. Reaction 
H4 is always important, whereas Reaction H5 is significant 
only at high ionization levels. 

Now Reaction H3 for H~ production can instead be 
rewritten as two different reaction chains: 

H3H4 2R + e'- 1 — > H2 e^ + 7 (A2) 
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H3H5 H + e" + H+ i — > 2H + 7 (A3) 

that occur at rates riYi^n^- R^Ri/ R^i- 
and nun^- R3R5 / Rii- , respectively. The first of these 
reaction chains forms molecular hydrogen using electrons as 
catalysts whereas the second leads to a reduction in the ion- 
ization level. 

We can treat in the same manner. Its decay channels 
are Reactions H7 and H8 which sum to a total destruction 
rate of 



R„+ ^nuRr + n^-Rs. 



(A4) 



Once again we have neglected photo-destruction which is 
important only at very high redshifts. 

Then Reaction 6 becomes the two reaction chains 



H6H7 
H6H8 



H 



2H-f 

- H+ + e" 



which are analogous to those for H 
plet starting with H2: 



Ha + H+ + 7 (A5) 
2H -f 7 (A6) 
There is a second cou- 



H9H7 
H9H8 



Ha + H+ + H 
H2 + H+-He- 



Ha H+ + H 
3H. 



(A7) 
(A8) 



Of these, H9H7 is the more important at low ionization levels 
and will strongly suppress destruction of Ha via collisions 
with protons. 

Using the reactions listed in Table |l], but replacing Re- 
actions H3 through H9 with the reaction chains derived 
above, we can write down an equation for the rate of change 
of Ha abundance: 



dnB_2 
dt 



"Ho 



R3R4, RaR? 
Rii- R^+ 



RgRs , r> I n 

n-H+^c- 1- nnRio + n^-Rii 



(A9) 



We can now write down the equation for the rate of 
change of Ha, then substitute for Jijj- and n-^+ . 



dnj- 



dt 



nHnii-R4 + nun„+R7 
2 

- riHa [nH+i?9 + nnRio + n^-Rii] 

a R3R4 , a ReRr , 
n-H ^o- 1" "-H — h nHnH2?iH+ 



(B4) 
R9R7 



Rtj 



= nn 



RaRa RcR7 
n^- — h 



Rv- 



(B5) 



- TlHo 



"h+'^c- -5 ^ "-Hrtio + n^-Rii 

Rti + 



Equation B5 is identical to Equation A9 derived above. 



APPENDIX B: DERIVATION USING 
EQUILIBRIUM VALUES OF H AND H+ 

The second derivation uses the fact that the destruction 
rates for H~ and H^ are high to derive equilibrium values 
for njj- and n„+ which can then be eliminated from the 
equations. 

Consider first the equation for the rate of change of 

= nnn^- R3 — n^- [nuRi + R5] . (Bl) 

The destruction rates within the square brackets are very 
large which means that rijj- will rapidly evolve to the equi- 
librium value in which creation and destruction of rijj- bal- 
- /dt X 0. Then we have 

nnn^-Rs _ nun^- R3 



n„+ 



nnRi -I- nH+i?5 -Rh- 
Similarly for n + we have 

riYiny^+Re + nH2?tH+-fi^9 

n^iRj + n^-Rs 
nnrij^+Re + nH2»H+-R9 



(B2) 



(B3) 



